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Abstract. We compute exactly the mean perimeter {L(T)) and the mean area 
{A{T)) of the convex hull of a random acceleration process of duration T in two 
dimensions. We use an exact mapping that relates, via Cauchy's formulae, the 
computation of the perimeter and the area of the convex hull of an arbitrary two 
dimensional stochastic process [x{t),y{t)] to the computation of the extreme value 
statistics of the associated one dimensional component process x(t). The latter can be 
computed exactly for the one dimensional random acceleration process even though 
the process in non-Markovian. Physically, our results are relevant in describing the 
average shape of a semi-flexible ideal polymer chain in two dimensions. 



1. Introduction 

Consider a set of N points witli positions {ri,r2, ■ ■ ■ on a two dimensional plane. 
How does one characterize the shape of this set? A standard characterization of the 
shape of this set is done by constructing the convex hull associated with this set. This 
is simply done in the following way. Let us imagine these points as nails stuck on a two 
dimensional board. Let us take a closed elastic or a rubber band, stretch it so that in 
includes all the points and then release it and let it shrink (see Fig. [T]). The elastic will 
shrink till it tightly encloses the points and can shrink no more (see Fig. Q. At this 
point, the shape of the elastic resembles a polygon and this is called the convex polygon 
or convex hull associated with the set. It is called convex because of the property that 
the line segment joining any two points on this polygon is fully contained within the 
polygon. 

Properties of convex hulls have been widely studied in mathematics, computer 
science and also in the physics of crystallography in connection with the so called 
Wulff construction. Convex hulls are widely used in computer aided image processing, 
in particular for pattern recognition [T]. Finding efficient algorithms to construct 
the convex hull of a set of points has led to numerous studies over the last few 
decades [21 El IH |5|, |6l [7] . In addition, convex hulls are used as a standard estimator for 
the home range for animal movements in ecology, leading to several important practical 
applications [H |9l [10]. For a recent review on the history and applications of convex 
hulls, see Ref. [12]. 
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Figure 1. Convex hull of a set of 9 points obtained by shrinking an elastic, initially 
stretched outside the set of points, up to the point where it can shrink no more. The 
resulting final polygon, shown by dotted red line, is the convex hull associated with 
this set of points. 

Of particular interest are the statistical properties of the convex hull associated 
with a set of random points on the plane. Consider, for instance, the set of N 
points whose coordinates {fi, f2, ■ ■ ■ j^n} are random variables drawn, in general, from 
a joint distribution P {ri,r2, ■ ■ ■ For each realization of the points, one can 

construct the associated convex hull and compute observables such as the perimeter 
L, area A or the number of vertices V of this convex hull. Clearly these observables 
change from one realization of points to another and are themselves random variables. 
Given the underlying distribution of the points P {fi,f2, ■ ■ ■ j^n), a challenging hard 
problem is to compute the statistics of the observables associated with the convex 
hull, e.g., the distributions such as P{L), P{A) or P{V). Even computing the mean 
perimeter (L) , mean area (A) and the mean number of vertices (V) on the convex 
hull is, in general, a hard problem. This is not just a mathematically challenging 
problem, but has several important practical applications, as reviewed recently in Ref. 
[12] . While several results are known (see [I2l for a review) in the case when the 
points are distributed independently and identically, i.e., when the joint distribution 
factorises, P (n, r2, • • • , "Hv) = Ili^=iP{ri), very few results are known when the vertices 
are correlated. 

A classic example when the points are correlated is the case when fj represents 
the position of a two-dimensional random walk at time step i. Each component of the 
position of the random walker evolves via the Markov evolution rule 

Xi = Xi-i + r]x{{) (1) 
Vi = Vi-i + T]y{i) (2) 
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starting from xq = yo = where the jump lengths rix{i) and riy{i) in the x and y directions 
are random variables, independent from each other and are also independent from one 
step to another and each is a zero mean Gaussian variable with a finite variance a^. The 
walk evolves upto steps and one can ask, for instance, what is the mean perimeter 
and the mean area of the convex hull associated with this random walk trajectory? 
These results are known in the continuous-time limit when N ^ oo and ^ keeping 
the product Na"^ = 2DT fixed, where T is the total duration of the walk and D is the 
diffusion constant. In the continuous-time limit, Eqs. ([l|-([2]) reduce to the Brownian 
motion in a plane 



starting at the origin x(0) = y{0) = 0, and where r]x{t) and riy{t) are Gaussian 
white noises with zero mean and the two-time correlators, {rix{t)ri,j.{t')) = 2D6{t — t'), 
{Vy{'t)Vy(t')) = 2D6(t — t') and {i]x(t)riy(t')) = 0. In this continuous-time limit, the mean 
perimeter of the convex hull associated to this Brownian motion of total duration T 
was first computed by Takacs, (L(T)) = y/lQn D T [13J. Note that while the \/DT 
dependence of the perimeter is expected due to the diffusive nature of the path, the 
computation of the prefactor a/IGtt is highly nontrivial. Later, the mean perimeter of 
the convex hull associated with a two dimensional Brownian bridge (where the walker 
is constrained to come back to the origin after time T) was also computed exactly by 
Goldman, (L(T)) bridge = Vn^DT [14]. The mean area of the convex hull for a free 
Brownian motion of duration T , {A{T)) = n DT, was first computed by El Bachir and 
Letac [151116], while the corresponding result for the bridge, (A(T))bridge = {2it/3)DT 
was computed only very recently [T7] . Once again, while the linear D T dependence is 
expected (as the area A ~ L^), the prefactors were nontrivial to compute. Moreover, 
the mean perimeter {Ln{T)) = anVDT and the mean area {An{T)) = l3n{DT) of the 
convex hull of n independent two-dimensional Brownian motions each of duration T was 
exactly computed recently for all n [TTl [12] both for free Brownian paths and Brownian 
bridges and the prefactors a„ and were found to have nontrivial n dependence. 

While the positions of a two-dimensional Brownian motion are correlated, the 
time evolution of the process is still Markovian. An important question is if one can 
compute the statistics of the convex hull associated to a two-dimensional non-Markovian 
stochastic process. This would then be a nontrvial generalisation. The purpose of this 
paper is to present exact analytical results for the mean perimeter and the mean area 
of the convex hull associated with the trajectory of a two-dimensional non-Markovian 
stochastic process, namely, the two dimensional random acceleration process. In this 
process, the position {x,y) of a particle in the plane evolves in continuous-time t via 



dx 

dy 
dt 



(3) 



(4) 



d'^x 
~d£^ 
d^ 

dt^ 



Vy(t) 



Vx{t) 



(5) 



(6) 
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Figure 2. Convex hull of 30 points generated using the discrete random acceleration 



process defined in Eqs. (30l-(31). The convex hull is constructed using the Graham 
span algorithm [5]. 



where rix{t) and riy{t) are zero mean Gaussian white noises as in Eqs. Q-Q, with 
two-point correlators, {r]^{t)r]^{t')) = 2D6{t - t'), {r]y{t)'r]y{t')) = 2D6{t - t') and 
{Vx{t)Vy{'t')) = 0- Henceforth, for simphcity, we will set D = 1 without any loss of 
generality. Note that due to the presence of the second derivative in Eqs. ([5])-([6]), the 
evolution is non-Markovian. If however, one defines the velocities, = dx/dt and 
Vy = dy/dt, then in the position- velocity phase space, the evolution equation involves 
only first derivative in time and hence the joint (r, v) process becomes Markovian, 

^^..W (7) 
^-,,W. (8) 

At time instant t the particle is identified by its position f{t) = {x{t),y{t)) and 
its velocity v(t) = {v.j.{t) = dx/dt, Vy{t) = dy/dt). We assume that the particle starts at 
the origin {x = 0,y = 0) at time with initial zero velocities: {vx{0) = 0,Vy{0) = 0). 
We evolve the process up to time T and then construct the convex hull associated with 
the trajectory in the (x — y) plane and compute the statistics of observables such as 
the perimeter L{T), the area A{T) and the number of vertices V{T). A particular 
realization of this process and the associated convex hull is shown in Fig. (|2|. 

Our main results are summarized as follows. We show that the mean perimeter and 
the mean area of the convex hull associated to a two-dimensional random acceleration 
process of duration T, starting from the initial zero velocity conditions Vx{0) = Vy{0) = 0, 
are given by the following exact expressions 




(L(r)) = \ —T'^'"^ = (2.1708 . . .) T^/^ (9) 
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= (0.100199. . .)Tl 



(10) 



As in case of Brownian motion, while the T dependence of the mean perimeter and the 
mean area easily follow from dimensional grounds, computation of the prefactors turns 
out to be rather nontrivial and they are the main results of this paper. We could not 
compute the mean number of vertices analytically, for which we present only numerical 
results. 

The random acceleration process is perhaps the simplest non-Markovian process 
and hence, no wonder, that it has been intensively studied both in the Physics and 
the Mathematics literature. Even though it is simple, the first-passage properties of 
this process are highly nontrivial even in one dimension, though several exact results 



are known [IHlliniEniiaiESlEniiaiESlETlEHlEgll^ There are 



several applications of this process, notably in the continuum description of semiflexible 
polymer chain [22], of fluctuating linear interfaces with dynamic exponent z = 4 [2S] and 
also in the description of the statistical properties of Burgers equation with Brownian 
initial velocity [31]. Thus, apart from being mathematically interesting, our results in 
this paper concerning the statistical properties of the convex hull of a two-dimensional 
random acceleration process are physically relevant in describing the average 'shape' 
of a ideal semi-flexible polymer chain in two dimensions (without the excluded volume 
interaction and with one end fixed at the origin while the other end is free). 

In deriving our exact results, we map the problem of computing the perimeter 
and the area of the two dimensional random acceleration process to the problem of 
computing the statistics of extremes of a one dimensional random acceleration process, 
following the general method introduced recently in [TTl [I2] . We remark that a similar 
mapping between the convex hull of a set of independent points drawn from a unit disc 
and an effective one dimensional random acceleration process was also noticed recently 
in the context of the so called Sylvester's question [35]. The Sylvester's question asks: 
given a set of independent points (drawn from a unit disc), what is the probabablity 
that all N points lie on the convex hull? In other words, in our notation the probability 
P{V = N) where V is the number of vertices on the convex hull. Based on the mapping 
to the one dimensional random acceleration process, the authors of Ref. [32] derived 
exact asymptotic results for P{V = N) for large N. 

The rest of the paper is organized as follows. In Section 2, we show how one can 
use Cauchy's formulae to map the problem of computing the perimeter and the area 
of the convex hull of a generic two dimensional stochastic process to the problem of 
computing the moments of the maximum and the time at which the maximum occurs 
for the associated one dimensional component stochastic process. We then focus on the 
random acceleration process and show how to compute explicitly the mean perimeter 
and the mean area of its convex hull using this mapping. In Section 3, we present the 
results of numerical simulations. We conclude in Section 4 with a summary and open 
problems. The detailed calculation of the moments of the maximum of a one dimensional 
random accelertion process, required for our main results, is relegated to the Appendix 
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A. 



2. Mapping to the extreme value statistics of a one dimensional process 
using Cauchy's formulae 

In Refs. \T2\ it was shown that the problem of computing the perimeter and the area 
of the convex hull of any two dimensional stochastic process [x{t), y{t)] can be mapped to 
computing the statistics of the value of the maximum and the time of occurrence of the 
maximum of the one dimensional component process x{t). This was achieved in pTj IT2] 
by exploiting the formulae for the perimeter and the area of any arbitrary closed, convex 
curve in two dimensions, derived originally by Cauchy [35] • We start this section by 
briefly reviewing this useful mapping and then use this mapping to compute explicitly 
the mean perimeter and the mean area of the convex hull of the two dimensional random 
acceleration process deflned in Eqs. Q-Q. 

Consider any closed convex curve C on a plane (see Fig. |3]). We choose the 
coordinate system such that the origin is inside the curve C. A key quantity, introduced 
by Cauchy [36], is the so called support function. Consider any direction 6. For flxed 
6, consider a stick perpendicular to this direction and imagine bringing the stick from 
inflnity and stop when it flrst touches the curve C (see Fig. |3|. At this point, the 
distance M{6) of the stick from the origin is called the support function in the direction 
6. It is clearly a very intuitive quantity as it measures how close can one get to the curve 
C in the direction 6 coming from inflnity. Mathematically speaking, consider all points 
on the curve C with coordinates [X(s),F(s)] (where s is the distance measured along 
the curve that just parametrizes the curve), consider their projections in the direction 
6 and then compute the maximum projection, i.e., 

M{e) = max [X{s) cos(0) + Y{s) sin(0)] . (11) 

Once this support function M{6) is known, then the perimeter L of C and the area A 
enclosed by C are given by Cauchy's formulae 



L = M{e)d9 (12) 







2 Jo 



M\9) - {M'{9)y\ de (13) 
where M'{9) = dM/dO. For example, for a circle of radius R = r, M{6) = r and one 



recovers the standard formulae: L = 2iTr and A = vrr^. The formulae in Eqs. (12)-(13) 
are generalisations of the circle formulae, valid for arbitrary closed convex curve C. A 
simple proof of the formulae can be found in the Appendix A of Ref . [T2j • 

Now, consider any two dimensional stochastic process [x(t),y{t)] up to time T, 
starting from [x(0) = 0,y{0) = 0]. For every realization, i.e., for every trajectory of 
the process, let us construct the corresponding convex hull C, which clearly varies from 
one trajectory to another. The idea then is to apply the Cauchy's formulae to each of 
these closed convex curves C and then average over all realizations to compute the mean 
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Figure 3. Support function M{d) in the direction 6 for a closed convex curve C on a 
plane. 



perimeter and the mean area. To apply Cauchy's formulae to a particular convex hull 
C, we need to first compute its associated support function M{9) defined in Eq. (11). 
It looks rather complicated to first compute the locus [X(s),F(s)] of the convex hull 
and then maximize over s. However, the crucial point to realize pTl IT2] is that M{6) 
represents the maximum projection of the convex hull C in the direction 6, which is also 
precisely the maximum projection of the trajectory [x{t),y{t)] itself in the direction 6, 
i.e., 

M{e) = max [X{s) cos{e) + Y{s) sin(^)] 

see 

= max^ cos(^) + y{t) sin(^)] . (14) 

So, in principle, if we have the list of cordinates visited by the trajectory up to time T, 
we can compute the support function M{6) by maximizing over time t. Let us denote 
Z0{t) = x{t)cos{6) + y{t)sm{6). For a fixed 6, zg{t) is just a one dimensional process 
parametrized by 6. Thus the support function M{6) is just the maximum of this one 
dimensional process ze{t) over the time interval t G [0,T]. Furthermore, assuming that 
the maximum of Z0{t) (with < t < T) occurs at time t^, we have 

M{e) = max [ze{t)] = ze{tm) = x(t„) cos(^) + y{tm) sin(^). (15) 

where tm implicitly depends on 6. The derivative M'{6) needed in the area formula in 
Eq. (13) can then be computed by deriving the right hand side (r.h.s) of Eq. (15) with 
respect to 9. This gives 

M'i9) = -x(t^) sin(^) + y{tj cos(0) + ^ ^ ^ ^ . (16) 

au at t—tm 
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Figure 4. Schematic picture of m{T), tm and y{tm)- 



However, since zg{t) is maximum at t = t^, by definition, dzg{t) / dt\t=t^_ 
the last term on the r.h.s vanishes and one simply gets 

M'{9) = -x{tm) sin(e) + y{t^) cos(^). 



and hence 



(17) 



The equations (15) and (17) are the essential ingredients behind the mapping of 
the original two dimensional convex hull problem to the statistics of the maximum of 
an effective one dimensional process. Now, an additional simplicity can arise if the 
two dimensional stochastic process is rotationally invariant, as in the case of Brownian 
motion and also for the random acceleration process. In this case, any average over the 
realizations is independent of the angle 9. Thus, the integral / = f^'" {F[M{9)]) d9 just 
reduces, thanks to the isotropy, to / = 27r(F[M(0)]), for any function F[M]. Setting 
9 



in Eqs. ([15j) and (17) we note that 

M(0) =x(t„ 
M'(0) = y{t^ 



m[T) 



(18) 
(19) 



where m{T) denotes the maximum of the one dimensional component process x{t) over 
the interval < t < T and t^ denotes the time at which x{t) achieves its maximum. 
On the other hand M'(0) is the value of the other independent one dimensional process 
y{t) at the time t = tm when the first process x{t) attains its maximum (see Fig. (|4])). 

Thus, the mean perimeter and the mean area, taking average over Cauchy's 
formulae, using isotropy, and Eqs. (jlsj) and (19), are simply given by 



(L(r)) = 2-K (M(0)) = 27r(m(T)) 
{A{T)) = vr Um\0)) - ((M'(0))^) 



TT 



{m\T)) - {y\t„^)) 



(20) 

(21) 



So far, the results in Eqs. (20) and (21) are quite general and hold for any arbitrary 
two dimensional rotationally invariant stochastic process [x{t),y{t)]. To compute the 
mean perimeter and the mean area, we then need to compute three quantities for the 
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associated one dimensional component process x{t), namely the first two moments of 
the maximum m{T) of the one dimensional process x{t) over the interval < t < T: 
(i) {m{T)) (ii) {m^iT)) and also (iii) {y^itm)) where tm is the time at which the process 
x{t) attains its maximum value. For the Brownian motion, these three quantities can 
be easily calculated [171 [12] . Here we show that for the random acceleration process as 
well, one can compute these three quantities explicitly. 

For a one dimensional random acceleration process a;(t), evolving via d'^x/dt^ = 
rix{t), over the interval < t < T, the Laplace transform of the distribution of the 
maximum was obtained recently by Burkhardt [32]. This Laplace transform has a 
rather complicated form (expressed as an integral over Airy function) (see Eq. (A. 4) in 
Appendix A). Computing explicitly the moments of m(T) from this Laplace transform 
turns out to be rather nontrivial. However, fortunately this can be done as we show in 
Appendix A. In particular, for the first two moments, we get 

(m(T)) = W — T^/' (22) 




(23) 



Finally, it rests to compute {y^itm))- This can be done in two steps. First, we note from 
the evolution equation, d'^y/dt^ = rjy{t) [where (fjyit)) = and {'Uy{t)r]y{t')) = 26(t — t')], 
that at any fixed time t, {y'^{t)) = (2/3) t^ for all t. This is easily obtained by integrating 
the evolution equation and computing {y'^it)) using the delta correlation of the noise. 
Secondly, averaging over the distribution p{tm\T) of the time of the maximum tm of x[t) 
one gets 



{y'itrr.)) = 3(4) 



dtm t^p{tm\T). 



(24) 



So, we need to know the distribution p{tm\T) of the time tm at which a random 
acceleration process achieves its maximum over the interval t G [0,T]. Fortunately, 
this is exactly known from our recent calculation 

;i-c). 



pitm\T)=C6iT-tm) + 



V2 



-3/4 



(T-t, 



1-1/4 



where the constant C has the exact value 



C = 1 




0.387628 . . . 



Using this result in Eq. (|24|) we finally get 

15 



{y\tm)) 




(25) 



(26) 



(27) 



Substituting the results from Eqs. (|22|), ([23|) and (|27j) in Eqs. (|20|) and ([2T|) 
immediately gives our main results 



rp3/2 ^ (2.1708. . .)T3/2 



(28) 
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3. Numerical simulations 

To verify our main theoretical predictions, we have also performed simulations of a 
discrete version of the random acceleration process in two dimensions. A particle at 
discrete time, n = 0, 1, . . . , A^, is identified by its position (x(n), y{n)) and its velocity 
{vx{ri), Vy{n)). At time n + 1 the velocity is updated by the discrete-time version of the 
Markov evolution rules in Eqs. ([T]) and ([S]): 

Vx{n) = Vx{n - 1) + r]x{n) 

Vy{n) = Vy{n - I) + riy{n) (30) 

where rj^ (n) , rjy (n) are uncorrelated Gaussian normal variables each with zero mean and 
variance o"^ = 2. The position (x(n), y{n)) is then updated: 



10 

(29) 



x{n) = x{n — 1) + Vx{n) 

yin) = y{n- 1) + Vy{n) (31) 

The process is started from x(0) = ?/(0) = fa;(0) = fj,(0) = and iterated for A^ — 1 
steps and then the convex hull is constructed using the Graham scan algorithm [2] . An 
example is shown in Fig. ([2]). The output of each realization is the ordered set of the 
V vertices sorted clockwise from {xi^yi) to {xv-,yv)- The perimeter and the area of the 
convex hull can then be computed using the formulae: 

y I 

L = ^\l{xi+i-Xi)'^ + {yi+i-yiy (32) 

1 ^ 

^ = -2 {^iVi+i - Xi+iVi) (33) 

i=l 

where we assume xy+i = xi and yv+i = yi- It is important to note that the minus 
symbol in Eq. (33) is due to the clockwise sorting of the points: indeed, the area 
formula (commonly called the surveyor's formula) gives the algebrical value of the area 
of a polygon; it gives the real area if the points are sorted counterclockwise, which is 
not the case here. 

The perimeter L, the area A and the number of vertices V are stochastic variables 
which fluctuate from one realization to another. We computed (L), (A), {V) for different 
values of A^. For each A^, we average over 10^ realizations. For large A^, the leading 
asymptotic behaviors of the mean perimeter and the mean area should be given by the 
continuous-time limit results in Eqs. (28) and (29) 

{L)=^N'/' + 0{N) (34) 



{A) = — 
^ ' 192 



^A^3^0(Ar5/2). (35) 
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Figure 5. Behavior of (L). Symbols: numerical data for the model defined by Eqs. 
(30) and (31 ). Solid line: continuum model prediction, ■\/37r/2 ^ 2.1708038 . . .. 
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Figure 6. Behavior of (A). Symbols: numerical data for the model defined by Eqs. 
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In Figs. (|5| and (|6|, we present our data respectively for the average perimeter (L) and 
the average area (A) of the convex hull, as a function of A^. Both the scaling and the 
numerical prefactor of the leading behavior in N are recovered. The agreement between 
the simulation and the analytical results is excellent. 

Furthermore, something rather interesting can be remarked on Fig. [7} For large 
A^, the leading behavior of the average number of vertices seems to be (V) ~ \/N. 
Numerically we find 

{V) r^CiVN + 02 + 0(1/ ^{N)). (36) 

A best fit gives Ci = 1.735 ± 0.01 and C2 = 2.6 ± 1.0. We were not able to derive this 
result analytically and it remains an interesting open question. 
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4. Summary and Conclusion 

In summary, in this paper we have computed analytically the mean perimeter (L(T)) 
and the mean area {A{T)) of the convex hull of a two dimensional random acceleration 
process of duration T, starting at the origin with zero initial velocity. We have used an 
exact mapping, via Cauchy's formulae, that relates the computation of the perimeter 
and the area of the convex hull of a two dimensional stochastic process to computing 
the extremal statistics associated with the corresponding one dimensional component 
stochastic process. Physically our results describe the average shape of a semi-flexible 
ideal polymer chain in two dimensions. 

We have also presented numerical results that suggest that the mean number of 
vertices V of the convex hull associated to the discrete-time version of the random 
acceleration process up to steps scales as {V{N)) ^ ci\fN + C2 + 0{l/\fN) for 
large A^. Proving this result and the computation of the constants Ci and C2 remain a 
challenging open problem. 

It is interesting to compare our exact results for the convex hull associated with a 
two-dimensional random acceleration process with those for the convex hull assocaited 
to a planar Brownian motion. In particular one notes that the ratio r = {A) / {L)"^ 
has very different values for the two processes. For a random acceleration process, 
using our exact results in Eqs. (9) and (10), we get r = (5/192)^2/3 = 0.02126.... 



The corresponding ratio for a free Brownian motion, using results mentioned in the 
introduction, is r = 1/16 = 0.0625 which is larger than the random acceleration process. 
The geometrical implication of this result is that the convex hull associated with the 
Brownian motion is, roughly speaking, more space filling than the one associated with 
the random acceleration process. In the latter case, the convex hull tends to be more 
'elongated' as evident from a typical representation in Fig. ([2]). This is somewhat 
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expected from the statistical mechanical consideration of the semi-flexible polymer 
chains described by the random acceleration process: the typical configuration of such 
a curvature-driven chain tends to be elongated as it costs more energy to bend or curve 
them as compared to a typical Brownian trajectory of a Rouse chain. 

There are several directions in which our work can be possibly extended. In this 
paper, we managed to compute only the mean perimeter and the mean area of the 
convex hull of the random acceleration process. It remains a challenging problem to 
compute the fluctuations or even the full distribution of these observables associated 
with the convex hull of the random acceleration process. The mapping we discussed in 
this paper is quite general and holds for arbitrary two dimensional stochastic process. 
It would be interesting to exploit this mapping to compute the mean perimeter and 
the mean area of the convex hull of other interesting two dimensional non-Markovian 
stochastic processes. For instance, for the process d"'f/dt"' = r/(t) (where n = 1 and 
n = 2 correspond respectively to the Brownian motion and the random acceleration 
process) with n > 2, it would be interesting to see if one can compute the statistics of 
the associated convex hulls. 

Appendix A. The moments of the mciximum value distribution 

Consider a one dimensional random acceleration process x{t) in the time interval [0,T]. 
At time t = the process starts from xq = with zero velocity, Vo = 0. Let 
m(T) denote the global maximum value achieved by x{t) during the interval [0,T]. 
We wish to compute the cumulative probability distribution of the maximum, i.e., 
Prob [m{T) < z] . Clearly, this is the probabihty that the value of the process, starting 
initially at x = 0, does not exceed the level z. In other words, Prob [m{T) < z] — 
Prob[x(t) < z, for all < t < T]. We then introduce the change of variable x{t) = 
z — x{t). The process x{t) is also a random acceleration process, starting from z (with 
zero velocity). The probability that m{T) is smaller than z can then be expressed in 
terms of the survival probability of the new process 5(t), i.e., the probability g(2;,0,T) 
that the process x{t) stays positive over the initerval t G [0, T] given that it started from 
the initial value 5(0) = z. We then write 

Prob [m(T) < z] = Prob[{a;(t)} <z\ = Vmh[{x{t)} > 0] = g(^, 0, T).(A.l) 

Note that g(0, 0, T) = and q{z oo, 0, T) = 1. 

The moments of the maximum m can be expressed in terms of the survival 
probability: 

/•oo rjn rOO fj POO 

(m"(r))= / dzz''^ = - dzz''^[l-q]^n dz z""-^ [1 - q(z,0,T)] (A.2) 
Jo dz Jo dz Jo 

where we have used integration by parts and the boundary condition q{z — )■ oo, 0, T) = 1. 

Let q{z, 0, s) denote the Laplace transform of the survival porbability: 

roo 

q(z,0,s)= dTe-'^ q(z,0,T). (A.3) 

^0 
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This quantity can be expressed as an integral [221 



- - q{z,0,s) = 



dF 



-^^Ai 



1 , V 2' 3 F 

40F 
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(A.4) 



where F |, is the incomplete Gamma function. We can thus introduce the Laplace 



transform of the moments in Eq. (A. 2) and obtain 



oo 



dTe-'^ {m''{T)) =n dz z 

Jo 

Performing the integral over z we get 

dTe-'^ {m''{T)) = r(n + 1) 



n— 1 



~ dF 



F5/3 



r 



1 2 s3/2 
2' 3 F 



4^/7^ 



POO 



dF 



Ai 



^2/3 



40F 



2s^ 
3 F 



Next we introduce a change of variable y = 2s^/^/(3F), 

r(n + 1) /3\"+i 



dTe-'^ {m'^iT)) 



S2 



|n+l 



dyy'^ 3 Ai 



f) 



2/3 



1 , 



40f 



As a result, the s dependence of the r.h.s is very simple, just proportional to s '^"/^ ^ 
and hence its inverse Laplace transform (with respect to s) can be done trivially to give 



(m"(T)) 



2 



r(|n) io 



dy 3 Ai 



I) 



2/3 



1 + 



r 



4^/7^ 



T2".(A.5) 



At this point, it turns out to be useful to express the Airy function in terms of the 
modified Bessel function K^{z) using the following identity [37] : 



Ai 



2/3' 



'IT 



3\t 



Then we have 



where 



(m"(T)) 



T(n) 



2j y37rr(|n) 



y'/'K^Ay). 



j{n) rp'^r 



dyy-KUy) 



1 + 



4^AF 



(A.6) 



(A.7) 



(A.8) 



We now show that the integral J*^"-' can be computed for any non-negative integer 
n. Starting from the definition of the incomplete F function, r{a,y) = dt t°'~^^ , 
we first use the following decomposition formula valid for any non-negative integer n. 
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which can be easily obtained by repeated integration by parts. We then spht the main 
integral J^"-* into three parts 1^"''^ = Ii + I2 + I3 with 



/•oo 

h= dyy-KM 
Jo 

h = i-irnn +h K. {y)T{-n - 1 

Z Jo Zti i Z 

k=i 2 Jo 2-K 



Appendix A. 0.1. Evaluation 0//2 Using the identity 

r(n + -)r(-n--,y) = — — / dt 



2' ' 2""' y"+V2yo t + y 

and the formula [371|, 

dy e-yK^{y) 7re*Ki(t) 



'0 y + t a/^ cos(z/7r) ' 

I2 can be recast: 

/, = (-1)- r dte-*r+i/2 r J^±LK^{y) 
^ Wo io 271^ t + y 3^^^ 

and the integral over y can be performed 

j2 = (-i)M rftrKi(t) = (-i)"Ji 

Jo 3 

Appendix A. 1. Evaluation of Ii 
The integral Ji can be computed [37] 



3/ v^r(^n) 



Appendix A. 1.1. Evaluation of Using the integral 

v/^r(/i + i/3)r(/i-i/3) 



JO 3 



2/^ r(/i + i/2) 

with n = n — k + 1/2 and the change of variable j = n — k we get 

-^3 = /— 5'n 



(A.IO) 



where 

^ _ r(j + 5/6)r(j + i/6)r(n - j + 1) 

This sum can be performed using Mathematica but it gives a rather cumbersome 



expression. Hence we prefer to keep this as a discrete sum in Eq. (A. 11). 
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Appendix A. 2. The final result 



Putting together all these integrals in Eq. (A. 7) we get 

" -1)" + 1 T{n) 



{m^T)) 



n - 
32 



1 



r(|n) , 



3n 

T^ 



(A.12) 



where Sn is given in Eq. (A.ll). For example, for the first four moments are given 
explicitly by 



(A.13) 
(A.14) 

(A.15) 
(A.16) 
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